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ABSTRACT 

We discuss a simple model for the growth of supermassive black holes (BHs) at the 
center of spheroidal stellar systems. In particular, we assess the hypotheses that (1) 
star formation in spheroids and BH fueling are proportional to one another, and (2) the 
BH accretion luminosity stays near the Eddington limit during luminous quasar phases. 
With the aid of this simple model, we are able to interpret many properties of the QSO 
luminosity function, including the puzzling steep decline of the characteristic luminosity 
from redshift z ~ 2 to z = 0: indeed the residual star formation in spheroidal systems is 
today limited to a small number of bulges, characterized by stellar velocity dispersions 
a factor of 2 — 3 smaller those of the elliptical galaxies hosting QSOs at z > 2. A simple 
consequence of our hypotheses is that the redshift evolution of the QSO emissivity 
and of the star formation history in spheroids should be roughly parallel. We find 
this result to be broadly consistent with our knowledge of the evolution of both the 
global star formation rate, and of the evolution of the QSO emissivity, but we identify 
interesting discrepancies at both low and high redshifts, to which we offer tentative 
solutions. Finally, our hypotheses allow us to present a robust method to derive the 
duty cycle of QSO activity, based on the observed QSO luminosity function, and on 
the present-day relation between the masses of supermassive BHs and those of their 
spheroidal host stellar systems. The duty cycle is found to be substantially less than 
unity, with characteristic values in the range 3 - 6 x 1CT 3 , and we compute that the 
average bolometric radiative efficiency is e ~ 0.07. Finally, we find that the growth in 
mass of individual black holes at high redshift (z^,2) can be dominated by mergers, and 
is therefore not necessarily limited by accretion. 

Subject headings: quasars: general - galaxies: nuclei - galaxies: active - black hole 
physics - accretion 
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1. Introduction 

The discovery of remarkable correlations between the masses of supermassive BHs hosted 
at the centers of galaxies and the global properties of the parent galaxies themselves (see, e.g., 
Magorrian et al. 1998; Ferrarese & Merritt 2000; Gebhardt et al. 2000; Graham et al. 2001) begs 
for interpretation. Several groups have noted the natural link between the cosmological evolution of 
QSOs and the formation history of galaxies (see, e.g. Monaco et al. 2000; Kauffmann & Haehnelt 
2001; Granato et al. 2001; Ciotti & van Albada 2001, Cavaliere & Vittorini 2002; Menci et al. 2003, 
and references therein) . The investigation of these interesting correlations looks promising not only 
to yield a better understanding of how and when galaxies formed, but also to obtain information 
about the QSO population itself (Ciotti, Haiman, k, Ostriker 2001; Yu & Tremaine 2002). For 
example, it may help us understand the well known but puzzling fact that the characteristic QSO 
luminosity (obtained from the QSO luminosity function, see, e.g., Pei 1995; Madau, Haardt & Rees 
1999, Wyithe & Loeb 2002) drops from z ~ 2.5 to z ~ by a factor of 35 ± 15. On the face of it, 
this result is surprising, since BHs can only grow, due to accretion or to mergers, and more massive 
BHs are expected to be more luminous on average, provided a sufficient amount of fuel is available. 

Here we focus on a few specific points raised by the general remarks above: (1) What drives 
the evolution of the steep decline with cosmic time of the quasar luminosity density, and of the 
characteristic quasar luminosity? (2) What is the expected relation between the cosmological 
evolution of the total emissivity in star-forming galaxies and that of the total emissivity of the 
quasar population? (3) How can one use scaling relations between the BH mass (hereafter Mbh) 
and the host galaxy properties to determine the QSO duty cycle at redshift z ~ 0? 

The rest of this paper is organized as follows. In § 2, we state our hypotheses and we list the 
observational inputs required by our approach. In § 3, we illustrate the technique adopted, and we 
explore quantitatively its consequences by linking the star formation history to the QSO evolution 
and applying it to explain the decrease of QSO mean luminosity with decreasing redshift. Then, 
in § 4, we present robust estimates of the QSO duty cycle and derive the mean accretion efficiency. 
Finally, in § 5, we conclude by summarizing the main results and the implications of this work. 

Throughout this paper, we adopt the background cosmological model as determined by the 
Wilkinson Microwave Anisotropy Probe ( WMAP; Bennett et al. 2003) experiment. This model has 
zero spatial curvature, and is dominated by cold dark matter (CDM) and a cosmological constant 
(A), with Q m = 0.29, ^ = 0.047, and 17a = 0.71, a Hubble constant Hq = 72 km s _1 , an rms 
mass fluctuation within a sphere of radius 8 h^ 1 Mpc of erg = 0.9, and power-law index n = 0.99 
for the power spectrum of density fluctuations (Spergel et al. 2003). These values are consistent 
with their determinations by most other methods (Bridle et al. 2003; Bahcall et al. 1999). 
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2. Basic Assumptions and Model Ingredients 

A widely accepted consequence of the so-called Magorrian relation, i.e., the (present-day) 
approximately linear relation between Mbh and Ms, the host spheroid stellar mass, is that the 
bulk of BH fueling in AGNs must be associated with star formation in the spheroidal components 
of their host galaxies (Monaco et al. 2000; Page et al. 2001; Granato et al. 2001; 2002; Cavaliere 
& Vittorini 2002). In this paper, we examine the simplest possible form of this association, namely 
the hypothesis that spheroid star formation and BH fueling are - at any time and in any system - 
proportional to one another with the proportionality constant independent of time and system. 

Since most of the mass of BHs appears to have assembled within a narrow redshift interval 
Az ss 1 around z « 2 (Boyle et al. 2000; Stoughton et al. 2002), in practice this hypothesis 
needs to hold only during this redshift interval, in order explain the local linear relation between 
BH and spheroid mass. One could argue that the energetic output from the forming central BH 
is the driving physical process that at the end will establish the galaxy mass (with the required 
proportionality). Alternatively, stars would form first, and then the BH is formed from reprocessed 
gas. In this case, a source of fuel for the BH growth with the required proportionality (namely 
mass losses from the newly formed stars), is available in a natural way. One can imagine that 
both of the above scenarios lead to a linear BH vs. spheroid mass relation at z = 0, but the 
strict proportionality in mass accretion rates into BHs and spheroids may not hold at all redshifts. 
Nevertheless, it is interesting to ask whether the simple hypothesis above is consistent with other 
observational data at both lower and higher redshifts, where some mass is still being added to both 
BHs ad spheroids, since this test can reveal information about the physical process of the BH and 
spheroid mass assembly. 

Since BH fueling should inevitably lead to some form of QSO activity, we make a second simple 
hypothesis, namely that the BH accretion luminosity always stays near the Eddington limit when the 
QSO is in the luminous, or "on" phase". This is apparently different from other proposals in the 
literature that are variants of the "feast or famine" model (Small & Blandford 1992), and which 
posit that QSO activity declines towards redshift z = owing, at least in part, to a significant 
decrease in the fueling rate (Cavaliere et al. 2000; Haiman & Menou 2000; Kauffmann &: Haehnelt 
2000). In fact overall "activity", i.e. the luminosity density evolution of QSOs, is the product of 
their characteristic number density Nq*(z) and their characteristic luminosity Lq*(z). The product 
Nq*Lq* may decline due to a decline in fueling that shuts off AGN activity and primarily leads 
to a decline in Aq*. But it is a separate question to ask what causes the surprising, but well 
observed decline with increasing time (Pei 1995; Boyle et al. 2000; Stoughton et al. 2002) in the 
characteristic luminosity Lq*(z). 

Coupled with the Magorrian relation, the above two hypotheses allow us to make several simple 
predictions, that will be described in detail in the following sections. Before we present our results, 
we list in detail the observational inputs required by our approach. 

The first observational input of our analysis is the present-day luminosity function (hereafter 
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LF) of spheroids, where the number of spheroids per unit volume with luminosities in the interval 
(L$,Ls + dL$) is defined to be given by $s(Ls)<ZLs- A composite LF was presented recently by 
Salucci et al. (1999), who considered the LF of four different types (E, SO, Sa/Sab, and Sbc/Scd) 
of galaxies separately, and inferred the total spheroid LF by assuming that on average, the spheroid 
components contribute 90%, 65%, 40%, and 10% of the light of the above galaxies, respectively. 
The composite spheroid LF is therefore represented by the sum of four different "Schechter-law" 
distributions 

&s(Ls) = J x ( 7 ) exp --— , (1) 

where log ($ S «/Gpc~ 3 ) = 5.89,5.95,6.03,6.45, log (L s «/L©) = 10.18,10.02,10.10,9.90, and a, L = 
0.95,0.95,1.0,1.3, for E and the bulges of SO, Sa/Sab, Sbc/Scd galaxies, respectively. Benson, 
Frenk, &: Sharpies (2002) have recently derived the spheroid luminosity function for a small sample 
of 90 bright field galaxies by decomposing the bulge and disk components, while Bernardi et al. 
(2002) and Sheth (2003) have computed (see also Yu & Tremaine 2002) the velocity function of 
early-type galaxies in the SDSS. While neither of these can serve as a substitute for the full spheroid 
luminosity function to replace equation (1) above, this should undoubtedly be possible in the near 
future by decomposing a large sample of fainter late-type SDSS galaxies into their bulge and disk 
components. 

The second ingredient is the quasar LF and its evolution with redshift, 

s (T x _ $Q*/Lq*{z) . (9 , 

Wq, z) - [Lq/Lq ^ )]A + [Lq/Lq.C*)]/* • {2) 

the optical data in the rest-frame B band can be well fitted by pure luminosity evolution, with the 
characteristic luminosity Lq* evolving with redshift as 



i e<*(l 



e 



L Q *{z) = L Q *(0)(1 + z)^" 1 J z + e £ Zm - (3) 

We adopt the fitting parameters given by Madau, Haardt & Rees (1999), Pi = 1.64, (3^ = 3.52, = 
1.9, £ = 2.58, £ = 3.16 and ckq = 0.5. The characteristic space density and luminosity are provided 
by Pei (1995) in a standard CDM cosmology with Hq = 50 km s _1 as log ($Q*/Gpc -3 ) = 2.95 and 
log [Lq*(O)/L ] = 13.03: we adopt these values with appropriate redshift-dependent re-scalings 
to our ACDM cosmology. 

Finally, the third ingredient is the Faber- Jackson relation (1976) 



0.62 
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1 300km s- 1 ) ' (4) 
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in the relatively more recent version of Davies et al. (1983), coupled with the Mrh — & relation 
(Ferrarese k Merritt 2000, Gebhardt et al. 2000, Tremaine and Yu 2002), 

10 9 M Q ~ V 300km s' 1 J ' U 
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Equation (4) is only approximately true, and the slope turns shallower for galaxies with velocity 
dispersions below 3 a < 170 kms -1 . Likewise, the exponent in equation (5) is currently under debate 
(Ferrarese & Merritt 2000; Gebhardt et al. 2000); here we refrain from a critical assessment of the 
different values found in the literature, and accept the slope ~ 4 as approximately the true value 
for both relations. Thus, to a good (but not necessarily perfect) accuracy, both the Mbh — c and 
the Faber-Jackson relations indicate a proportionality to the fourth power of the central velocity 
dispersion, implying the following linear relation: 

~ 0.016—. (6) 

When expressing equation (6) above in term of galaxy mass instead of luminosity, it is found that 
the implied median BH mass fraction to stellar mass is 0.13% of the mass of the bulge (Kormendy 
& Gebhardt 2001). 



3. Linking the Star Formation History and the Evolution of Quasars 

As emphasized above, it is widely believed that the bulk of BH fueling in AGNs must be 
associated with star formation in the spheroidal components of their host galaxies. In this section, 
we examine the hypothesis stated in § 2 above, namely that spheroid star formation and BH fueling 
are - at any time and in any system - proportional to one another with the proportionality constant 
independent of time and system. Under the assumption that quasars radiate a fixed fraction e of 
their accreted mass, an obvious consequence is that the redshift evolution of the QSO emissivity 
and of the star formation history in spheroids should be roughly parallel to each another. As we 
shall see, we find this result to be broadly consistent with our knowledge of the evolution of both the 
global star formation rate, and of the evolution of the QSO emissivity, but we identify interesting 
discrepancies at both low and high redshifts, to which we offer tentative solutions. 

The evolution of the total UV luminosity density in stars at 1500 A (galaxy rest frame) with 
redshift is given (in a standard CDM cosmology with Hq = 50 km s -1 ; Madau & Pozzetti 2000) by 



This is related to the total star formation rate density as 

PS 



^' UV Moyr^Mpc" 3 (8) 



10 27 



for a Salpeter IMF (Madau et al. 1998). Figure 1 shows (dashed curve) this star formation 
rate density (SFRD), with and appropriate redshift-dependent re-scaling to our adopted ACDM 



3 For example, a fit to galaxies in Virgo cluster by Dressier et al. (1987) gives a mean exponent of ~ 3.5, while 
Bernardi et al. (2003) found instead a value of ~ 4 for the exponent. 
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cosmology. This SFRD is close to that derived more directly in the recent work by Porciani & 
Madau (2001). 

The evolution of the total rest-frame B-band luminosity density in quasars can be obtained 
from equations (2) and (3) as 

/•oo 

jq,b= / dLL$ Q (L,t). (9) 
J o 

Under the assumption that quasars radiate a fixed fraction e of their accreted mass (see discussion 
below), this is related to the total BH accretion rate density (pq,b, hereafter BARD) as 

A ho \ec 2 ( p QyB \ 
jQ ~lM 0r iMpc-3j' (10) 

where ^4boi = H-2 is the bolometric correction, and e = 0.071 is the radiative efficiency, derived 
in the next section. Figure 1 shows (solid curve; displaced upward by a constant factor of 770 for 
clarity) the BH mass accretion rate density, re-scaled to our adopted ACDM cosmology. 

As is well known, both the SFRD and the BARD exhibit a steep rise from z = 0to2; = l — 2, 
a peak at z ~ 1 — 2, and a decline towards still higher redshifts. This is broadly consistent with 
their expected parallel evolution under our simple set of assumptions. Both the SFRD and the 
BARD still have significant observational uncertainties. While the steep decline at low redshift 
is relatively secure, the current SFRD and the BARD determinations could both turn out to be 
underestimates at high redshifts, due to yet-undetected populations of galaxies or AGNs (e.g. due 
to dust obscuration). A critical review of the uncertainties is beyond the scope of this paper; we 
here simply take the current determinations at face value, and examine discrepancies at both low 
and high redshifts from our simple model. 



3.1. Low Redshifts (z < 2) 

3.1.1. Why Does the Characteristic QSO Luminosity Evolve? 

We start by demonstrating that the bulk of BH formation, and consequently the bulk of QSO 
activity, must have occurred in galactic systems dominated by massive, luminous bulges. In fact, 
the spheroid light distribution is known to approximately satisfy a Schechter-like distribution (see, 
e.g., eq. 1) 

/(^ocf^V'expf--^ (11) 



.Ls*J V ^s*, 

with a « 1.2 ± 0.1 (Salucci et al. 1999; Benson et al. 2002; Bernardi et al. 2003). Then, from 
Mbh oc Lg , it follows that one half of the mass in BHs is in systems with luminosity Lg / Lg* > £1/2, 
where £1/2 is defined by 

r a+1 exp(-£)d£ = K (12) 
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and for a = 1.2, this yields i\j2 — 0.5. In a more detailed computation, using the composite 
spheroid luminosity function given in equation (1), we find that the luminosity above which half 
of the integrated light is emitted corresponds to a spheroid with luminosity Mb ~ —20.5, only a 
factor of ~ 2.5 fainter than the luminosity of the well-known giant elliptical M87 (Mb = —21.42). 
Thus the bulk of the mass density of BHs reside in giant ellipticals, and, if the current situation is 
not anomalous, the bulk of the growth of SMBHs must also have occurred there ( or in progenitor 
systems). We also know that the bulk of star formation in spheroidal systems took place as early as 
redshift z > 2, as indicated, for example, by the mean stellar ages in ellipticals (Hogg et al. 2002; 
Bernardi et al. 2003), and of bulge populations (e.g. Proctor et al. 2000; Ellis et al. 2001), and 
the Butcher-Oemler or Gunn-Dressler effects (Margoniner et al. 2001). 

At present, the disks of spiral galaxies dominate the global star formation rate (Fukugita et 
al. 1998; Benson et al. 2002; Hogg et al. 2002), and the mean age of stars in spiral systems is 
perhaps a factor of two younger than that in spheroidal systems. It follows, given our hypothesis, 
that BH growth in the local universe is dominated by relatively small bulges that live in galaxies 
denominated as spirals. Fortunately, the hypothesized relation between nuclear activity and star 
formation can be directly tested at low redshift. For example, Percival et al. (2001) obtained 
morphological information for the host galaxies of nine bright, (My < —25.5) QSOs, classifying 
six of them as "disks" , and the remaining three as "spheroids" . The bulk of the local population 
of identified QSOs live in disk dominated systems. The sample studied by Percival et al. was 
approximately three magnitudes brighter than the characteristic luminosity of the local population 
of QSOs (Pei 1995), and almost all low luminosity AGNs are known to reside in spiral systems. 
Our conclusion would therefore likely be strengthened by a survey going to magnitudes fainter than 
studied in the Percival et al. work, closer to Lq*. 

It is natural to ask what the consequences would be of the hypothesis that BH fueling, when 
it does happen, stays near the Eddington limit. This assumption is not unrealistic: in fact, for 
a handful of nearby AGNs, the BH masses can be directly estimated by reverberation mapping 
(or by a cruder "photoionization method"; see Wandel et al. 1999), and for these sources, the 
Eddington ratio can be directly inferred. From the 19 nearby AGNs listed in Table 3 and in Figure 
5 of Wandel, Peterson, & Malkan (1999), one derives L = (0.01 — 0.3)LEdd f° r the Seyfert 1 objects, 
while the two QSOs have L/LEdd = 0.2 and 0.3. The bolometric correction for very hard and 
IR/submillimeter radiation, using the mean quasar spectrum as given by Elvis et al. (1994), is 
about a factor of ~ 3, which would bring the luminosities of the two QSOs in the Wandel et al. 
sample close to Eddington limit. 

The assumption of always maintaining the Eddington luminosity, if applied to an individual 
BH, predicts an increasing QSO luminosity Lq (in the "on" state), due to the trivial fact that for 
every BH the mass is monotonically increasing. However, this is not necessarily in conflict with 
observations that describe the evolution of the characteristic luminosity for a population of quasars. 
Clearly, if all galaxies remained equally active, then the mean luminosity inferred from the observed 
QSO luminosity function would increase, but if the typical member of the population changes with 
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time the naive expectations may be incorrect. 

Let us simply assume as an empirically verified fact that at the present day, star formation is 
active primarily in disk dominated systems. A decrease of a factor 20-50 in the characteristic quasar 
luminosity Lq* would then be naturally obtained by combining the Faber-Jackson (L oc a 4 ) and 
Magorrian (Mbh oc f 4 ) relations with a reduction in the characteristic central velocity dispersion 
(a) in the hosts of QSOs by a factor in the range 2.1-2.7. A decrease of this amount is quite 
natural, when one considers the mean (luminosity weighted) central velocity dispersion associated 
with the ellipticals at redshift z = 2 (a ss 400 km s -1 ), and that associated with spiral bulges at 
redshift z = (a ss 200 km s _1 ), as derived by the Faber-Jackson relation. This argument thus 
provides a straightforward interpretation of why the typical QSO luminosity decreases from z = 2 
to z = 0. Furthermore, there is explicit observational evidence (Thomas et al. 2002) that large 
ellipticals are older than small bulges of spirals, supporting the decrease in a towards z = as the 
reason behind the decrease in the characteristic quasar luminosity. 

A slightly less steep drop in the central velocity dispersion could be acceptable, by simulta- 
neously allowing for quasar luminosities to decrease to somewhat sub-Eddington values towards 
low redshifts. The latter scenario is consistent with the Eddington ratios of the two quasars in 
the Wandel et al. sample. Applying the reverberation mapping technique to an extended quasar 
sample (e.g. selected from the Sloan Digital Sky Survey, SDSS) would distinguish directly among 
these two options. 

3.1.2. Why Does the Total Quasar Emissivity Evolve? 

We next consider whether the observed steep evolution of the total quasar emissivity (or 
BARD) is consistent with the star formation rate density (SFRD). Figure 1 shows the BARD and 
SFRD. However, the SFRD includes contributions from both disk stars and spheroid stars: here 
we discuss corrections to this diagram to obtain the SFRD in spheroids alone. 

Corrections for disk star formation are likely to become large at low redshifts. The fraction of 
the total stellar luminosity density at z = contributed by stars in disks vs. stars in spheroids (the 
latter including both the bulges of spirals and ellipticals) has been estimated by several authors. 
Fukugita et al. (1998) and Hogg et al. (2002) find that spheroids contribute ~ 40% of the luminosity 
density (but Benson et al. 2002 find spheroids to contribute significantly less than this fraction). 
Furthermore, it is well-known that the stellar populations in present-day spheroids are old, and 
most models place their formation epochs at z > 2. Nevertheless, we are interested in the amount 
of ongoing starformation in these spheroids at z = 0. The lower limits on the ages of the spheroid 
stellar populations come from various methods; one of these is the colors. B — V and V — I 
magnitudes can be determined to an accuracy of ~ 0.1 mag (e.g. Ellis et al. 2001), and these colors 
typically change by ~ 1 mag in a Gyr of evolution (Leitherer et al. 1999). It follows that < 10% of 
the present-day luminosity in these systems can arise from young stars. In turn, this implies that 
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less than ~ 5% of the total SFRD at z = is occuring in spheroid systems (i.e. in the bulges of 
spirals). On the other hand, at redshifts of z = 2 — 4, the Lyman break galaxies are believed to be 
large bulge systems in the process of formation, and the fraction of the starformation seen at these 
redshifts associated with bulges essentially unity. 

The Butcher-Oemler (also known as Dressler-Gunn) effect can also be used to derive the 
contribution of spheroids to the total SFR as a function of redshift in galaxy clusters. The BO 
effect then gives the rate of expected decline in the emission from SMBHs, since these tend to be 
in high mass ellipticals which are well represented in the BO clusters. Observational work on the 
BO effect has shown that the number fraction of blue galaxies in clusters, fo, increases with z in a 
linear relation (see, e.g., Newberry et al. 1988, Andreon & Ettori 1999, Metveier et al 2000). For 
example, in one of the most detailed works, based on the analysis of 295 POSS-II clusters with 
redshift < z < 0.4 (Margoniner et al. 2001), the authors found 

/ 6 ^1.3 ±a5 z + c, (13) 

where c is a small additive constant, of the order of 0.02 ± 0.01 There are important caveats in 
using the BO effect for the spheroid correction. First, there are large variations in from cluster 
to cluster (e.g. the z = 0.83 cluster studied by Van Dokkum et al. (2000) has an estimated 
fb = 0.22 ± 0.09). It is also not clear at the present time where the star formation responsible for 
this blue light is occurring. While this blue light may represent ongoing star formation in spheroids, 
Abraham et al. (1996) argue that the blue light arises in the disks that flare up as spirals fall into the 
cluster potential. However, interpreting the blue fraction as the fraction of ellipticals undergoing 
starformation, the BO effect would support the general conclusion that spheroids contribute only 
a few percent of the total starformation rate in the present-day universe, while this fraction rises 
steeply towards higher redshifts. 

There are other promising methods to estimate the spheroid contribution to the total SFRD. 
For example, one could measure local starburst activity in the dense obscured centers of galaxies 
in the infrared bands, and identify this with the local star formation rate in bulges. It should 
also be possible to measure an accurate age distribution of stellar populations in the bulges of 
spirals, as well as in elliptical galaxies, in large samples of SDSS galaxies, and hence to directly 
infer the time dependence of the star formation rate in spheroids as a function of redshift. An 
accurate measurement of the age distribution has already been achieved for a red subsample of 
SDSS galaxies (Jimenez et al. 2003). 

In summary, we here estimate the rough fraction of the observed starformation rate that is 
associated with spheroids at each redshift by multiplying the total SFRD shown in Figure 1 by a 
factor / sph = 0.05 + 0.95(z/2.2) 2 at z < 2.2. This ensures a smooth transition in the total SFRD 
being dominated by ellipticals at high redshift to it being dominated by disks of spirals at z = 0, 
with only residual starformation in the bulges of spirals, consistent with the arguments above. In 
Figure 2, we show (long-dashed curve) the corrected SFRD. As the figure demonstrates, including 
this correction improves the fit, in the sense of making the SFRD resemble the BARD more closely. 
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However, intriguingly, it does appear that the decline in the spheroid formation rate from z = 2 
to z = is too large, by a factor of ~three, when compared to the decline in the BARD. If this 
discrepancy holds up in future data, it would imply that the nuclear black holes can be fueled long 
after the star formation in the bulge has ceased. Since the bulge star formation has likely used up 
all the gas initially present in the bulge, the fuel would have to arrive from elsewhere. 

A simple assertion is that old stars that had formed in the bulge keep returning a fraction of 
their mass in winds. These stellar winds may provide dense, shocked material that can dissipate 
and serve as a fuel for the central BH. The mass-loss rate in winds in a starburst evolves as oc i _L3 
(where t is the time elapsed from the burst; see, e.g., Ciotti et al. 1991; Leitherer et al. 1999). We 
here use the wind mass-loss rate M w j n d = 1-5 x 10 -11 Lb£[" 5 M Q yr _1 between the ages of 0.5 and 
15 Gyr for a one solar mass model starburst galaxy, where t\s is the time elapsed from the burst in 
units of 15 Gyr, and Lb ~ 0.03 is the B-band luminosity of the model galaxy (in units of L ) at 15 
Gyr (Bruzual & Chariot 2000). Under these assumptions, ~ 80% of the stellar mass is eventually 
returned to the ambient medium. 

For our purposes, we regard mass lost in winds as new material available to fuel the central BH. 
It is clear that the total mass return rate from a passively evolving stellar population cannot accrete 
onto the central BH (otherwise BH masses will be two or three order of magnitude larger than those 
observed); nor can it all turn into stars (star formation in present day ellipticals is not detected 
at the level that would be implied). In order to solve these problems Ciotti & Ostriker (2001) 
performed numerical simulations of radiative feedback modulated accretion flows onto a SMBH at 
the center of a " cooling flow" galaxy. They showed that only a few percent ( or less ) of the available 
gas lost in winds effectively accretes onto the central BH, while the accretion luminosity during 
short episodes of bursts stays near the Eddington value (a similar conclusion would follow in the 
case of mechanical feedback; e.g. Tabor & Binney 1993; Binney & Tabor 1995; Binney 1999). 

Here we assume that a fraction 1.3 x 10~ 3 of the mass in winds accretes onto the central BH; i.e. 
the same fraction we had assumed for the "original" infalling gas earlier in this paper. This fits in 
well explicitly with the fractions inferred in Ciotti & Ostriker (2001). In Figure 2, we show (dotted 
curve) the total mass loss rate generated in winds from spheroid stars as a function of redshift. The 
thick solid curve (turning into the dot-dashed curve at high redshift, see discussion below) shows 
the total BARD inferred from the SFRD after mass loss from winds are added. We conclude that, 
if a significant fraction of the wind material ends up fueling the central BH, this brings the BARD 
and SFRD into quite reasonable agreement (to within a factor of two at all redshifts). 

We emphasize that our treatment in this subsection is phenomenological, and complementary 
to theoretical semi-analytical models (Haiman & Menou 2000) of the cosmological evolution of the 
QSO luminosity, or models based on Monte Carlo realizations of dark matter "merger trees" (see, 
e.g., Kauffmann & Haehnelt 2001). These models have found that to reproduce the observed decline 
in the QSO luminosity density, the "efficiency factor" for the fraction of gas accreted by the BH in 
a merger must decline towards z = 0. Our proposal here is radically different: instead of "starving" 
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the BHs in each galaxy, the QSO luminosity density drops due to the empirically-inferred drop in 
the formation rate of spheroids, and their bias towards smaller systems at lower redshifts. 

3.2. High Redshifts (z > 2) 

As at z < 2, at redshifts exceeding the peak of quasar activity (z > 2), the evolutions of the 
SFRD and BARD are, in fact, not parallel (see Fig. 1). It must be noted that the observational 
determinations of both quantities are much less certain at these high redshift than at low redshifts. 
The presence of a population of high-redshift, dust obscured quasars could, for example, reconcile 
the SFRD and BARD curves in our simple model. There is already evidence for such a population 
that could significantly increase the inferred BARD (e.g. Fabian & Iwasawa 1999); there is also 
some evidence that, unlike the optical LF, the soft X-ray quasar luminosity function stays flat out 
to redshifts z ~ 4 (Miyaji et al. 2001). We here simply take the current determinations at face 
value, and examine physical reasons that would explain the apparent discrepancy, if it holds up in 
future data. 

3.2.1. What Steepens the Evolution of the High-z Quasar Emissivity? 

One possibility is that the fueling rate of quasars is suppressed by intrinsic physical limits to the 
rate of accretion. Models in which the BHs shine with approximately their Eddington luminosity 
can naturally explain the observed evolution of the QSO luminosity function (Haiman & Loeb 1997; 
Haehnelt, Natarajan & Rees 1998; Wyithe & Loeb 2002), by associating the rise from z = 6 to z = 2 
with the increase in the nonlinear mass-scale in hierarchical cosmologies. Ciotti & Ostriker (2003) 
suggested that at the high characteristic densities at z > 3, Bremsstrahlung opacity may effectively 
limit the mass accretion rate onto a BH to a small fraction of the usual Eddington value. This 
idea is attractive because it provides a physical reason for the suppression of the fueling rate, and 
because the additional opacity may be relevant only at high redshifts, allowing "normal" accretion 
at z < 3. In fact, the data reviewed in section § 3.1.1 are consistent with L max = O.lLEdd (and a 
modest correction for beaming), but the effects of correspondingly increasing the Eddington time 
by a factor of 10 are only important at high redshifts, where it then becomes comparable to the 
age of the universe. 

In Figure 2, we show (dot-dashed curve) the evolution of the emissivity for a BH, L oc M cx 
M oc exp(/Edd*/*Edd) under the assumption that the hole grows exponentially on a timescale 
which is /j^j = 1/0.07 = 14 times the Eddington time tEdd = e x 4.6 x 10 8 yr. A multiplicative 
constant (1O _3,1 M0 yr _1 Mpc~ 3 for the curve shown in Fig. 2) can be used to represent the summed 
emissivity or accretion rate density of all quasar BHs, all of which are assumed to grow at the same 
rate from z = oo. As the figure reveals, the suppression of the accretion rates in all BHs to 7% 
of the Eddington value would naturally result in the observed steep slope of the quasar emissivity 
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evolution between 3 < z < 6, while not preventing star-formation to occur in a more extended 
spheroid region around the black hole. While attractive, this explanation suffers from a drawback, 
namely the fact that if all BHs can accrete only at 10% of the Eddington rate at z > 3.5, then 
their e-folding time will be ~ 5 x 10 8 years, making it apparently difficult to explain how the large 
(few xlO 9 M ) BHs in the SDSS survey were built by z = 6, when age of the universe is 8 x 10 8 
years, less than twice the e-folding timescale (Haiman & Loeb 2001). As we shall see below, this is 
less of a problem than might be expected, since the growth of individual SMBHs at high redshift 
is dominated by mergers, and not by accretion. 



A different, but potentially important ingredient in determining the relative evolution of the 
SFRD and the BARD at high redshifts is the importance of mergers (see also Volonteri et al. 2003). 
We next demonstrate that at high redshifts, the build-up of the mass of an individual BH is likely 
dominated by mergers between BHs. Such mergers may not have any effect on the total quasar 
emissivity (one can imagine merging all BHs in pairs, resulting in new BHs twice as massive as 
the original set, while preserving the accretion rate per unit BH mass). However, if mergers are 
frequent, one can imagine that this may help explain the high redshift discrepancy between the 
SFRD and the BARD. For example, one can imagine that a merger event at high redshift delivers 
new gas and triggers star-formation (but may not be able to increase the accretion rate onto BHs, 
per unit BH mass, if this quantity is already Eddington limited). 

The growth rate of BHs due to merging can be obtained from the characteristic dark matter 
halo number density as follows: 



Here we define the nonlinear dark matter halo mass-scale M m at redshift z as a(M n \)g(z) = 1, 
where a(M) is the r.m.s. mass fluctuation in spheres of mass M, and g(z) is the growth function 
at redshift z. For simplicity, we define the space density of halos as A m = M m x dN/dM n \, where 
dN/dM is the usual (comoving) halo mass function, adopted here from Jenkins et al. (2001), 
evaluated at M m . Under this assumption, N Q \ = N n \(M(t),t), and we have the following time 
derivative: 



3.2.2. The Growth of an Individual Black Hole due to Mergers vs. Accretion 



M, 



merg — 




(14) 



N nl = M n i 



d?N 



+ M nl 



dN 



(15) 



dM ni dt 



The first term on the right hand side vanishes by definiton (d 2 N/dMdt ~ at the nonlinear 
mass-scale). As a result, we find that 
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Mmerg = 0.13 X 0.1 X 1.3 X 10~ 3 X M al — ~ 1.7 X 10~ 5 X M nl . (16) 

This last result reflects the fact that without any accretion, the individual BH masses would 
grow only by coalescence of the BHs during halo mergers, and therefore the typical BH mass would 
simply track the nonlinear dark halo mass-scale. In order to describe the growth of BHs, rather 
than that of halos, we have assumed in equation (16) that a fraction 0.13 of the total mass in 
each halo is baryonic, the mass of stars is 10% of the baryons, and the mass of the central BH is 
1.3 x 10~ 3 that of the stars. 

We next find the growth of an individual BH due to accretion, using the Madau & Pozzetti 

(2000) star formation rate density, as follows: 

M acc = 1.3 x 10" 3 x (17) 

where ps is the comoving star-formation rate density as given by equations (7) and (8) above, and 
we take 1.3 x 10 -3 for the ratio of the BH mass to the spheroid mass from Kormendy and Gebhardt 

(2001) . 

Figure 3 shows the mass growth rates due to merging (solid curve) and accretion (dashed curve) . 
From this figure, we learn that at redshifts (2 < z < 4), the growth is dominated by mergers, while 
at low redshift (z < 2), the growth is proportional to the star formation rate. 4 According to 
the figure, at very low redshifts (z < 0.5), mergers are again important; however this regime is 
unphysical because of the so-called "over-merging" problem in the Press-Schechter formalism: the 
nonlinear mass-scale grows to that corresponding to clusters of galaxies; the galaxies, however, 
may preserve their identities, and hence it is no longer clear that the growth of BHs by coalescence 
in merging galaxies tracks this mass-scale. 



4. Radiative Efficiency and Duty Cycle of AGNs 

In this section, we define and derive the radiative efficiency and the duty cycle of AGNs; 
quantities that served as inputs in the previous sections. In the interest of clarity, let us first 
consider a population of N g identical galaxies over the Hubble time tn, each of which today (i.e. at 
t = tu) hosts a spheroid of mass Ms, and a BH of mass Mbh- Let us further assume that during 
the entire time elapsed from to tn, each BH had only two states: it was either "on" or "off". We 
identify the "on" state as the active quasar phase, and we define the duty cycle /q as the fraction 
of the time each BH spends in the "on" state. At any given time, the number of active quasars is 



4 Note that at high redshifts (z > 4), the mass buildup by mergers actually exceeds ~ 10% of the Eddington rate, 
in accordance with § 3.2.1 above. 
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then Aq = /qA^. In the "on" state, the BH grows by accretion at the rate Mrh, and shines at the 
(bolometric) luminosity Lq with a radiative efficiency e, defined as the fraction of the rest mass 
energy of the infalling gas converted to radiation. The remaining fraction (1 — e) of the rest mass 
then leads to the growth of the BH mass (Yu & Tremaine 2002). A simple algebra shows that 

e _ fQtnL Q N g _ t H L Q N Q 



1-e M^ u c 2 M BH N g c 2 M BU N g c 2 ' 

Here c is the speed of light; the numerator represents the total light emitted by all BHs, and 
the denominator represents the total mass in BHs today. In the third equality, we have used 
-^Q = /Q-Wg- Note that the last term involves only quantities that are, in principle, directly 
observable, and that it is independent of the duty cycle (Soltan 1982). Equation (18) describes the 
entire galaxy population, but a similar equation applies to individual galaxies: Lq/q^ = eMbhc 2 . 
This last expression does have a dependence on the duty cycle, which can therefore be written as 

/Q -^-^I^ (19) 



The simple toy model above demonstrates that (i) the radiative efficiency can be obtained 
independently of the duty cycle, and (ii) that the duty cycle can be obtained two different ways, 
based either on the number or on the characteristic BH mass of quasars. While the former method 
is conceptually more straightforward, as we shall see below, the latter avoids the divergence in the 
number of quasars and galaxies (due to the steep observed slope of the luminosity functions at the 
faint end). 

The step forward to a more realistic situation is to allow a distribution of galaxy and BH 
masses, and corresponding quasar luminosities, and to allow the BH masses and luminosities to 
change with time. In this case, equation (18) can be straightforwardly generalized to obtain the 
global average radiative efficiency (Soltan 1982). The total energy output of all quasars over all 
times per unit volume is given by 



rZmax ^£ /*Lmax 

u Q = -A hol dz— dLL$ Q (L,z) (20) 

JO az JL min 

where $q(L, z) is the observed quasar luminosity function, which we take from equation (2) in the 
B-band, A^ \ is the bolometric correction, which, from the composite quasar spectrum in Elvis et 
al. (1994), we find to be A^ i = L tot /L,B = 11.2, and dt/dz is obtained from the time-redshift 
relation in our chosen ACDM cosmology. Note that the integral in equation (20) converges both 
from above and below in luminosity, and in redshift, so that it is insensitive to the integration limits 
£min, L max , and z max . We find 

^ = 1.9 x 10 4 M Mpc~ 3 (21) 
The total remnant BH mass in the present universe is given by 
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p BH = f dMM$ BU (M) = / BH f dMM<Z> s (M), (22) 
Jo Jo 



where <3?bh and <3?s are the present day BH and spheroid mass functions, and /bh is the BH to 
spheroid mass ratio. We here adopt the spheroid luminosity function from equation (1), and convert 
it to a spheroid mass function using the M — L ratio for each of the 4 different types of galaxies 
in Salucci et al. (1999). Assuming further a constant BH-spheroid mass ratio /bh = 1-3 x 10~ 3 
(Kormendy & Gebhardt 2001), we find 

Pbu = 2.49 x 10 5 M Mpc" 3 . (23) 

It is worth noting the contribution to this mass density from the (E, SO, Sa/Sab, and Sbc/Scd) 
galaxies separately, which are, respectively, 1.19, 0.63, 0.50, and 0.17 M Mpc~ 3 . In other words, 
the bulges of late type galaxies are not a negligible contribution to the total spheroid mass density 
(and they dominate the total mass density for low spheroid masses). It is also worth noting that 
the present-day total spheroid mass density we obtain by integrating the spheroid formation rate 
in Figure 2 (which was based on a corrected version of the evolution of the total star formation 
rate in Madau & Pozzetti 2000) is in good agreement (to within 30%) of the value we find from 
summing up the inferred local spheroid densities from Salucci et al. (1999; see eqn. 1). 

Finally, from the ratio of equations (21) and (23), we find e = 0.071. Our result is in good 
agreement with recent work by Yu & Tremaine (2002), who find an average efficiency of e = 0.077. 
It is interesting to note that this agreement holds separately for the QSO light and BH mass 
density. Yu & Tremaine use the recent 2dF quasar luminosity function of Boyle et al. (2000), and 
a bolometric correction of 11.8, to find uq/c 2 = 2.1 x lO 4 M0Mpc -3 . They combine the recent 
Mbh — & relation of Tremaine et al. (2002) with the SDSS velocity function of early type galaxies 
in Bernardi et al. (2002), and make corrections for the additional BH mass in spiral galaxies and 
for the scatter in the Mbh — cr relation, yielding /Obh = 2.5 x 10 5 M Q Mpc -3 . 

The duty cycle defined in equation (19) is less trivially generalized for a population of evolving 
galaxies. Nevertheless, if we assume that the duty cycle does not vary with time (but we allow it 
to be a function of luminosity), then we can still obtain the duty cycle explicitly by comparing the 
present-day space density of quasars and galaxies. An immediate complication is that both space 
densities (see eqns. 1 and 2) diverge at the faint end. We therefore proceed by defining the average 
duty cycle of all quasars above luminosity Lq(x, 0), 

</Q - N>i = ' { ' 

where Lq(x, 0) is such that QSOs at redshift z = brighter than Lq(x, 0) emit a fraction x of the 
total quasar light Lq = J °° dL L$q(L, 0), and likewise, Mbr{x) is such that all BHs more massive 
than Mbh(^) sum up to the same fraction x of the total BH mass at z = 0, i.e. to xM^ n . 
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We may similarly generalize the definition of the duty cycle based on the characteristic BH 
mass (cf. the last term in equation [19]), by applying it to an individual present-day BH, as follows: 

, f \ _ ec 2 M BH (x) 

{jQ,M/x = r oo T , — 7777 > [ 2b ) 

Jo L Q (x,t)dt 

In principle, the time-integral in the denominator on the right-hand side must be taken over 
the (typical) luminosity history of the present-day black hole of mass Mbh- If mergers play an 
important role in the growth of this BH, then the integral must be performed separately, and 
summed over each branch along the "merger tree". In practice, we do not know this merger 
history, and instead we simply define the time-dependent luminosity Lq(x,£) such that QSOs at 
cosmic time t brighter than L,Q(x,t) emit a fraction x of the total quasar light at that epoch, 
Lq = Jq 00 dL L$Q(L,i). This definition assumes only that a monotonic relation is maintained 
between Mbh and Lq at all times, and it gives the correct duty cycle in the limit that mergers do 
not dominate the mass growth of individual BHs (note that the bulk of the total quasar light is 
emitted in a narrow peak around redshift z = 2.2; ~ 50% of the time-integral in equation (25) is 
contributed between 1.6 < z < 2.8). As we showed in the previous section (see Figure 3), mergers 
are likely important at z > 3 but do not dominate the growth of individual BHs at lower redshifts. 

The duty cycles obtained from both methods are listed in Table 1. We find that (/q,n)o.i — 
0.008 and (/q,n)o.9 — 0.05. Most importantly, our two definitions above do not guarantee that 
equations (24) and (25) give the same values. Here we find that the two methods agree well on 
the high mass end, while (/q,m) is systematically lower by a factor of ~two towards the low-mass 
end. One interpretation of this finding is that the most massive BHs gain nearly all their mass by 
accretion, while mergers contribute a comparatively larger fraction of the mass of lower mass BHs. 

Our results for the duty cycle being at the percent level is in good agreement with theoretical 
expectations (Ciotti & Ostriker 1997, 2001, Yu & Tremaine 2002), it is also similar to the values 
derived in Haehnelt, Natarajan, & Rees (1998), who obtained a QSO lifetime of £q ~ 10 7 yr at a 
Hubble epoch of ~ 10 9 ~ 10 yr or, in terms of the duty cycle, /q = tQ/tmibbic — 10~ 2 — 10 -3 . Similar 
quasar lifetimes can be independently derived from the spatial clustering of quasars (Haiman & 
Hui 2001; Martini k Weinberg 2001). 

Before we conclude this section, we stress a well known puzzle presented by the difference in 
the QSO luminosity functions in the optical and X-ray bands, and its consequences for our analysis. 
We repeated the derivation of the accretion efficiency from the X-ray luminosity function (Miyaji et 
al. 2001), and we find e ~ 0.045, about two times lower than we obtain from the B-band (applying 
in both cases a bolometric correction from Elvis et al. 1994: L tot /Lx = 38.1, and L to t/L,B = 11.2). 
If all QSOs emitted intrinsically with the Elvis spectrum, the above two numbers should agree. 
However, they differ by a factor of two. What does this mean? The simplest resolution is to 
assume that QSOs emit a universal spectrum, but the mean value of their flux ratio Lx/Lb is a 
factor of 71/45 = 1.6 times smaller than for the Elvis et al. sample (which consists of 47 unobscured 
QSOs). This says Lx/Lb ~ 0.18 . However, under the assumption of a universal spectrum, the 



total number of quasars 



/•oo d y rco 

N(>L mm)=i * dz dL {L , Z) (26) 

should then be equal in the optical and in the X-rays, if one uses the appropriate lower limits, 
£min,x=0.18L m i nj B- We find that this leads instead to Nx = 2Mb- In other words, under the 
assumption of a universal spectrum, the flux ratio Lx/Lb can be derived two ways: either from 
the total number, or from the total light, of quasars. With the published LFs, these methods give 
Lx/Lb ~ 0.5 and Lx/Lb ~ 0.2, respectively. This proves that quasars as a population cannot 
have a universal spectrum, a result that one could have obtained directly by comparing the X-ray 
and optical LFs, which have different shapes and redshift evolutions. A possible resolution is that 
the duty-cycle in the X-rays is ~ 3 times longer than in the optical. This would be supported by 
the recent results of Barger et al. (2001), who found that a large fraction of optical galaxies are 
Chandra AGN sources, implying a long X-ray activity cycle, ~0.5 Gyr. 



5. Conclusions 

In this paper, we discussed a simple, empirically based model for the growth of supermassive 
black holes (BHs) at the center of spheroidal stellar systems. Motivated by accumulating evidence 
for the strong link between the formation of spheroids and BHs, we hypothesized the simplest possi- 
ble form of this connection, namely that star formation in spheroids and BH fueling are proportional 
to one another, at all cosmic epoch and in all spheroids, regardless of their size. 

The main conclusions that arise from this hypothesis (augmented with a few other reasonable 
assumptions) are as follows. This simple model accounts for the puzzling steep decline of the 
characteristic luminosity of quasars from redshift z ~ 2 to z = 0: the residual star formation in 
spheroidal systems is today limited to a small number of bulges, characterized by stellar velocity 
dispersions a factor of 2 — 3 smaller those of the elliptical galaxies hosting QSOs at z > 2. We 
explored a very simple consequence of our hypothesis: the redshift evolution of the QSO emissivity 
and of the star formation history in spheroids should be roughly parallel to each other. We find 
this result to be broadly consistent with the evolution of both the global star formation rate, and 
of the evolution of the QSO emissivity, both of which exhibit a peak at redshift z ~ 2. However, a 
closer look reveals possibly interesting discrepancies at both low and high redshifts. 

At low redshifts, the spheroid formation rate, obtained by making simple corrections to the 
total star formation rate, appears to decline by a factor that is ~ 3 times larger than the decline 
in QSO emissivity. A possible solution we note to resolve this discrepancy is fueling of quasar BHs 
at low redshifts by the mass lost in winds from a passively evolving stellar spheroid population, 
formed at earlier epochs. A tentative discrepancy also exists at high redshifts (z > 2, beyond the 
peak of QSO activity), where the evolution of the star formation rate appears significantly flatter 
than that of quasar emissivity. While a population of hitherto undetected, obscured AGN at high 
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redshift (with the obscured fraction increasing towards high z) could resolve this discrepancy, we 
offered an alternative, physical explanation: quasar fueling rates at high redshift are limited to a 
fraction ~ 10% of the Eddington accretion rate. This limit depends linearly on the characteristic 
BH mass, and would therefore imprint a steep evolution of the quasar luminosity function as the 
characteristic mass-scale builds up exponentially. We also note that the masses of individual black 
holes at high redshift are not limited by accretion (at the Eddington or some modified Eddington 
rate), since we find that mergers dominate over accretion in determining the growth of objects at 
epochs z>2. 

Given our demographic assumptions, we compute the average duty cycle - the fraction of time 
SMBHs spend in the on state - as (3 — 6) x 10~ 3 , depending on BH mass, and we also find the 
mean bolometric radiative efficiency, e = 0.071, when averaged for the entire SMBH population. 

The considerations on this paper are tentative, but empirically based. It should soon be 
possible to considerably tighten constraints on the simple picture of coeval formation of BHs and 
spheroids: the redshift-evolution of the starformation rate in spheroid systems should be derivable 
by a more detailed analysis of the SDSS galaxy sample out to at least z ~ 0.5. 

We thank Martin Haehnelt, Martin Rees, Alvio Renzini, Scott Tremaine and Qingjuan Yu 
for useful discussions. L.C. was supported by MURST, contract CoFin2000, and by grant ASI 
I/R/105/00. 
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Table 1. 



X 


M BH (x) 




L Q (x,0)/L s * 


N>(x,0) 


(/q,n}ie 


(/q,m)z 


0.1 


4.33e8 


3.92e-5 


1.51e0 


3.02e-7 


0.0077 


0.00349 


0.2 


2.67e8 


1.14e-4 


0.80e0 


1.05e-6 


0.0092 


0.00404 


0.3 


1.82e8 


2.29e-4 


0.49e0 


2.33e-6 


0.0102 


0.00452 


0.4 


1.27e8 


3.94e-4 


0.30e0 


4.40e-6 


0.0112 


0.00517 


0.5 


8.91e7 


6.28e-4 


0.18e0 


7.85e-6 


0.0125 


0.00616 


0.6 


6.03e7 


9.69e-4 


0.93e-l 


1.40e-6 


0.0145 


0.00784 


0.7 


3.76e7 


1.50e-3 


0.42e-l 


2.67e-6 


0.0179 


0.0109 


0.8 


2.00e7 


2.40e-3 


0.14e-l 


6.03e-5 


0.0251 


0.0180 


0.9 


7.05e7 


4.48e-3 


0.20e-2 


2.19e-4 


0.0488 


0.0433 



Note. — The columns show respectively, the BH mass, the BH space density, 
quasar luminosity, local QSO space density, and duty cycles (last two columns, 
computed using two different methods) , for galaxies that contain a fixed fraction 
x (from column 1) of quasar light and both spheroid and BH mass. 
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Fig. 1. — The redshift evolution of the total star formation rate density (SFRD; dashed curve) 
and the total black hole accretion rate density (BARD; solid curve). The SFRD was adopted from 
Madau & Pozzetti (2000); while the BARD is obtained from the optical quasar luminosity function, 
assuming a bolometric correction of A^ Q \ = 11.2 and a constant radiative efficiency of e = 0.071 
(independent of redshift and quasar luminosity). The BARD is displaced upward by a constant 
factor of 1/e = 770 for clarity of presentation. 
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Fig. 2. — Modifications to the SFRD and BARD in Figure 1 that would bring the two quantities 
to parallel each other; as required by our simple set of assumptions. The modifications include 
(1) a correction for the spheroid vs. total SFRD (long-dashed curve); (2) the possibility of fueling 
BHs via stellar winds (dotted curve); and (3) a suppression of the BARD at high redshifts due to 
an extra source of opacity (dot-dashed curve). The thick solid curve (turning into the dot-dashed 
curve at high redshift) shows the BARD inferred from the SFRD after these corrections are taken 
into account; it tracks the BARD inferred from the optical quasar luminosity function (thin solid 
curve) . 
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Fig. 3. — The growth rate of an individual black hole with the characteristic BH mass at each 
epoch, from accretion (dashed curve) and mergers (solid curve), as a function of redshift. Also 
shown is the growth rate corresponding to accretion that is limited at all times to a fraction 0.07 
of the Eddington rate for the characteristic BH mass at each epoch. 



